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Abstract 

This talk discusses the effective field theory view of deep inelastic scattering. In such an 
approach, the standard factorization formula of a hard coefficient multiplied by a parton 
distribution function arises from matching of QCD onto an effective field theory. The 
DGLAP equations can then be viewed as the standard renormalization group equations 
that determines the cut-off dependence of the non-local operator whose forward matrix 
element is the parton distribution function. As an example, the non-singlet quark split- 
ting functions is derived directly from the renormalization properties of the non-local 
operator itself. This approach, although discussed in the literature, does not appear to 
be well known to the larger high energy community. In this talk we give a pedagogical 
introduction to this subject. 



1 Introduction 



Electron-nucleon Deep Inelastic Scattering (DIS) has played an important role in our under- 
standing of the strong interactions. By now, it is a standard topic in particle physics and 
quantum field theory books [U [2] . 

The cross section for this process can be written in term of two structure functions. These 
obey a factorization formula (see for example [TJ [3]) 
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The parton distribution function (PDF) <pi/h is often referred to as the "probability to find 
a parton i in hadron h" . Such a description, although very intuitive, requires a field theoretic 
definition. What do we mean by a "distribution function"? 

This question was answered in the early 1980's [U [S]. The quark PDF is 
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where n light-like vector, n 2 = and W n is a light-like Wilson line, 
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and = A a ^t a , with t a in the fundamental representation. 

Let us discuss ([2]) in more detail. First, in the definition of the PDF, an average over color 
and spin is understood. Second, since ipi and ipi are not at the same point, we must include a 
"gauge- link" to render the operator gauge invariant. This is the reason that the Wilson lines 
appear in the definition of the PDF. Notice that the product of the two semi-infinite Wilson 
lines can be written as finite Wilson line connecting the point and tn. Third, the PDF is just 
the Fourier transform of the diagonal matrix element of a standard dimension-three operator. 
Since the operator is non-local, its Fourier transform is a function and not a number. Fourth, 
the anti-quark PDF and gluon PDF can be defined in a similar way to the quark PDF. For 
example, the gluon PDF is defined with the quark field replaced by gluon field strength and 
the Wilson line is in the adjoint representation. For concreteness in the following we will focus 
on the quark PDF. 

The PDF satisfies the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equation 
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where P^ are the so called "splitting functions". These are usually calculated via two main 
methods. 

The first is based on the interpretation of the splitting function as the "probability" of the 
quark to "split" into a quark and a gluon [6j. At 0(a s ) one finds 
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where Cp = 4/3 for QCD. The full expression is in fact 
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and the singular terms 
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defined in the appendix, and 5(1 — £) must be added "by 



hand". 

The second method, often called the "OPE" method [H [8], is based on the fact that 
moments of the PDF are matrix elements of local operators. Thus one defines the moments 
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which satisfy a local renormalization group equation (RGE) 

d<fi n , 
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The splitting function is related to 7„ via 



7« = - / dxx n ~ x P k 
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There is also another approach possible, which is less known, and can be described as the 
"effective field theory" method. Going back to the factorization formula for F 2 , we can write 
it schematically as 



F 2 = CMaO ® </>(/•* ) + power corrections 

t t (6) 

Wilson coef. Operator 

The factorization formula can be interpreted as a result of matching QCD onto an effective 
field theory. G 2 is the Wilson coefficient extracted in the process, and the PDF is the matrix 
element of the operator in the effective theory. The "power corrections" correspond to the 
contribution of power suppressed operators. This is the standard expression one finds in a 
generic effective field theory approach. The only complication in this case is the non-locality 
of the operator. The scale \i on which the PDF depends is then just the cut-off scale of the 
effective theory. It is not surprising then that the cut-off dependence of the PDF is determined 
by an RGE which is schematically 



dln/i 

This is none other than the famous DGLAP equation, where the splitting function is simply 
the anomalous dimension. Since the PDF is a function, the anomalous dimension is a function 
too, and not a constant. 
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If the splitting function is just the anomalous dimension of the PDF, one can calculate 
it as one usually does in field theory. In particular, if we use dimensional regularization and 
a mass independent scheme, the anomalous dimension can be extracted from the 1/e pole of 
the loop corrections to the PDF. The only unusual aspect of the calculation arises from the 
non-local structure of the PDF. We will discuss this calculation in detail in the next section. 

This calculation was first performed, to the best of our knowledge, by Braunschweig, Horejsi 
and Robaschik [9]. But their paper is almost unknown. By the end of 2009 it had 5 citations 
on SPIRES [10J. For comparison, by the end of 2009 reference [5] had 3933 citations while 
references [7] and [8] had 589 and 919 citations, respectively. The situation is slightly better 
then these citation counts indicate, but still the third approach is much less known, especially 
to the general high energy community. 

The effective field theory approach has received much interest in the recent years, since the 
effective field theory one is matching onto is the soft collinear effective theory (SCET) [TTJ [T21 
[13]. Calculations of the anomalous dimension of non-local operators are standard in SCET. 
The application of SCET to DIS, as well as other hard QCD processes, was first discussed in 
P3]. In that paper the matching coefficients were calculated at tree level. Although SCET 
is the appropriate EFT for DIS, the effect of the soft gluons completely cancel for generic 
x O [T5]. In this case SCET is just "boosted QCD" and we can calculate the anomalous 
dimension by using QCD Feynman rules instead of SCET Feynman rules. We will demonstrate 
this approach by calculating P qq for the non-singlet case. The method will be technically 
different from [9], but conceptually very similar. For another pedagogical discussion of the 
calculation of the evolution (RGE) equations for non-local operators see appendix G of [IS] , 
where the calculation techniques are demonstrated for both covariant and light-cone gauges. 

2 Example of a Non-Local Renormalization: 
Non-singlet Splitting Function 

2.1 Feynman rules 

We begin by deriving the Feynman rules. We will use the definition of the PDF in equation 
(T5]), and replace the hadron state by a free quark. Since we are interested in the cut-off (UV) 
dependence of the PDF which is independent of the IR physics, such a replacement is justified. 
The "zero-gluon" Feynman rule is 






n ■ p — n ■ k) u(k) — u(k). 
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Figure 1: Feynman diagrams 



Taking the external states to carry momentum p we find that the zeroth order expression is 
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Averaging over the spins we would find — 1), i.e. the quark is carrying all of the momentum 
of the hadron. The zero-gluon Feynman rule is therefore 
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S(£ n ■ p — n ■ k) — 



In order to find the one-gluon Feynman rule we need to expand the Wilson lines in (j3J). A 
simple calculation yields 




t a g \6(£ n ■ p — n ■ k) — 5(£ n-p — n-k + n-l) 
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2.2 Feynman integration 

In order to find P qq we need to calculate three diagrams, see Figure [Q To demonstrate the 
techniques, we will calculate the leftmost diagram in Figure [TJ We will use Feynman gauge in 
the calculation. 

First we need to define light-cone coordinates. We define two light-cone vectors n 2 = 
0, n ■ n = 2. In particular we will use n = (1,0,0,-1), n = (1,0,0,1). Notice that 
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n ■ p = Po + P3 > 0. Any four vector can be decomposed as 

or = n ■ a— + n ■ a— + a]_, 
which also defines the "±" coordinates. In these coordinates d d k can be written as 



d d k = - dn ■ k dn ■ k d d 2 kj 
2 



We can now write down the Feynman integral, 
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where p (k) is the external (internal) quark momentum and 

Num. = (n • kf | + n ■ k Jt± - k]_ -. 

The total diagram is given by 

D\ = -ig 2 C F ^ u(p) ri^u(p) (10) 

The only difference from the usual Feynman integral is the presence of the delta function. 
It also forces us to split the integration according to the light-cone coordinates. Our strategy 
will be to calculate the n ■ k integral using residue theorem, to calculate the k± integral in 
d — 2 dimensions, and use the delta function to perform the integral over n ■ k |17j . 

In order to use the residue theorem for the n ■ k integral we need to find the poles of the 
integrand. They are given by 

k 2 + iO = n ■ kn ■ k + k\ + iO = n-k = — — 

n ■ k 

(k-p) 2 + i0 = n-k = n-p+ v , J± . (11) 

n ■ k — n ■ p 

In order to get a non-zero result we must ensure that the two poles are on the opposite sides 
of the real n ■ k axis. Otherwise we can close the contour in the half plane that contains no 
poles and get zero. There are two options then, 

n ■ k 

I) n-k>0 U n-k-n-p<0 =3- 0<n-k<n-p =3- 0< <1 



II) n ■ k < U n ■ k — n ■ p > =^ n ■ p < n ■ k < 0. 



n-p (12) 



Recall, though, that n ■ p > 0, so the second option is ruled out. Combining the first option 
with the delta function <5(£ n • p — n ■ k) we find that < £ < 1. Notice also that we have 
recovered the usual support property of <t>i/h{C)i namely, that the momentum fraction carried 
by the quark must be between and 1. 

After performing the n ■ k integral using residues, we perform the integral over k± using 
dimensional regularization in d — 2 dimensions. Due to the numerator structure, we can split / 
into 3 integrals that differ by their Dirac structure. In order to find the anomalous dimension 
we only need the Dirac structure that appear in the definition of <pi/h{C)-> l - e i 1 - The other 
Dirac structures lead to finite terms and do not contribute to the anomalous dimension. 



^^We can verify this point using another argument. The integral can be organized by the so called method 
of regions [HI QjJ] , which is also related to our effective field theory interpretation. The method of regions 
allows us to divide the momentum integration to various regions in which the integral momentum scales in a 
particular way. The sum of all the regions is equal to the full integral. 

For our integrals there is only a collinear region, namely a region in which the components of the momentum 
k scale as 

n.*~0(l), n-k~0(A 2 QCB ), ^~0(A QCD ). 
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We will therefore ignore the other Dirac structures in the numerator. In order to regulate 
the IR divergences we keep p 2 7^ 0. Performing the n ■ k integral for this structure by using 
the first pole in (ITTj) . we find 



1 

~2 
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k 2 + i0 
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(k-p) 2 + i0 



5(£ n ■ p — n ■ k) 



d d - 2 k ± 1 
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where we have defined r = n ■ k/n ■ p. The limits of integration are determined by the n ■ k 
integral as explained above. Performing the relevant k± integration, we find that the divergent 
part of the diagram is 



div 



X 



•— %^ u(p) rtiAp) (V) £ (i - e)r(e)(V - ,-o)- 

dr (1 — r)[r(l — r)]~ e <5(£ — r), 



Expanding in e we finally find 

1 1 C F a 



div 



f/ ,. y , [7r «(p)|«(p)2(i-owe(i-e). 



(14) 



(15) 



The other two diagrams can be calculated using the same methods. Using the one-gluon 
Feynman rule we find that 



Do = D? 



— ^u(p) iu(p) (4^ r(e)(V - .0)- 
n ■ p Atx 2 

e~ e (i - iY^eii - 0^(0 - 5(1 - b(2 - 6, - e ) 



(16) 



In order to expand in e we need to use the identity [20] , which is also proven in the appendix, 
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Expanding in e we finally find that the divergent part is 

1 



D 



div 



div 



+ C(e). 
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We also have 



d 4 fc~0(A 4 QCD ), fc 2 ^0(A| CD ) ^I~A 
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Since the anomalous dimension can depends only on n ■ p, it must be an 0(1) quantity, and so does /. 
We therefore need the numerator to scale as Aq CD , so only the k\ integral will contribute to the anomalous 
dimension. This integral involves only the i/i Dirac structure. 
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Notice that the singular terms 



and 5(1 — £) arise naturally in the course of the 



calculation and are not added "by hand" . 

We also need the wave function renormalization constant Zq [21] in Feynman gauge 
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The total expression for the matrix element at one loop accuracy is 



(19) 



Atai = Df v + Df^ + Df v + Z Q D 
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Averaging over the spins we find that 

</w(0 = S(l - + ^ \ 9(0 9(1 - [ (1 + a + 1*0-- 0] + finite + («s) 

(21) 



2.3 Renormalization 

We now calculate the RGE equation in the MS scheme in the "standard way" , see for example 
We define the renormalization factor Z(x,0 by 



w (x) = / dtz(x,o r n -(o = z(x,o ® r n -(x), 



(22) 



and expand Z in a s /47r, such that Z\ n i denotes the coefficient of the (a s /47r) n term. At tree 
level we have 4>^ e (x) = (f) T ^(x) = 5(1 — x), so 



Z M (x,0 = ^-a;) = i«5(l-|). 

In order to calculate the one loop expression for the renormalization factor we use 

<$f (x) = Z [0] (x,0 ® fffiO + Z {1] (x,0 ® 0[ o c f-(O 

to find 
or 

= - f ) + " f )'(f ) 7 1 p i? (f ) + 



(23) 



(24) 
(25) 
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where P qq is defined in (JSJ). In order to calculate the anomalous dimension we now follow the 
standard procedure [22] , treating Z and 7 as infinite dimensional matrices. Define 

O bare = Z ® O rcn - and — = 7 ®O ren -, (26) 

a ln/i 

and expand Z as 

00 1 

2 = 1 + ^!^. (27) 



k =i e 



The anomalous dimension 7 is given by 



7 = 2a s - — . 28 



In this way we find that 



7(x, = - 7>G 7^(7)+ (29) 



and the RGE is finally 



d</>(x) f 1 f 1 d£ 
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(/ln// !i ^ 1 (x,om = J jP&($<I>\t)+0(<4), (so) 

where we have changed £ — >• x/£. We have recovered the standard DGLAP equation at 0(a s ). 



3 Conclusions 

We have argued that the standard features of DIS, namely, factorization into a short distance 
coefficient times a parton distribution function and DGLAP evolution of the PDF, can be 
understood in the language of an effective field theory. In this interpretation QCD is matched 
onto an effective field theory. The short distance coefficient is simply the Wilson coefficient and 
the PDF is the matrix element of the operator in the effective field theory. The only difference 
from more familiar effective field theories is that the operator is non-local and therefore its 
matrix element is a function and not a number. The cut-off dependence of the matrix element 
can be determined in the standard way. In particular, by using dimensional regularization and 
a mass independent scheme, the anomalous dimension can be extracted from the 1/e pole of 
the loop corrections to the PDF. The resulting RGE equations are none other than the famous 
DGLAP equations. This approach, although discussed in the literature, does not appear to 
be well known to the larger high energy community. 

We have demonstrated this approach explicitly by calculating the one-loop correction to 
the non-singlet quark PDF. We have then derived the evolution equation from the divergent 
piece of the one-loop correction. This procedure can be generalized in several ways. First, we 
can calculate higher order corrections in a s . Second, we can include operator mixing, such as 
the calculation of P qg and P gq . Third, we can include power corrections. For example, one can 
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calculate in this method the anomalous dimension of twist 4 operators, which is still an open 
problem. 
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A The plus distribution 

The "plus" distribution — is defined such that for any test function fix) 

(1 -x)+ 



Notice that 



(31) 



^ f{x) 

x) ■ 

More generally [^(x)] is defined as 

J\xf{x) [h(x)] + = j\x (f{x) - /(l)) h(x). (32) 
When multiplying + by a function g(x) which is regular at x = 1 we have 

1 dxf(x)g(x) [h(x)} + = j\x [f{x)g{x) - f{l)g{l))h{x). (33) 



[g(x)h(x)] + = g(x) [h(xj\ + -S(1- x) J 1 dx [g{x) - g{\)) h(x). (34) 
Another useful property is the identity 

(l- x )»-J—- = (l-x)»-i. (35) 

[1 - X) + 

Thus we can write 

pS = ( 1 '^ + 2 ^jt^o; + 1 6{1 ~ = {1+e) jr^T + l 5{1 ~° = (t=j ) + - (36) 

In order to prove the identity 

6(x)6{l - x){\ - x)- 1 -" = - x) + (jz^j + 0(e), (37) 

we multiply (1 — x)~ 1 ~ e by a test function, f(x) and integrate 

dxf(x)(l-x)- 1 -^ [ dx[f{x)-f{l)]{l-x)- 1 ^ + f{\) [ dxil-x)- 1 -* (38) 

Jo Jo 
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Both terms on the right hand side are regular at x = 1, so we can expand in e and find 

1 dx f{x) {1 - x)- 1 -' = f dx fiyX ] - -101 + 0(e). (39) 
'o Jo 1 x e 

Notice that 

9(x)9(l - x) g(x) (1 - x)- 1 - = -i 5(1 - x) g(l) + g(x) (j^j + O(e). (40) 
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